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We analyze the thermal fluctuations of particles that have a short-range dipolar attraction and 
a long-range repulsion. In an inhomogeneous particle density region, or "soft phase," filamentary 
patterns appear which are destroyed only at very high temperatures. The filaments act as a fluctuat- 
ing template for correlated percolation in which low-energy excitations can move through the stable 
pattern by local rearrangements. At intermediate temperatures, dynamically averaged checkerboard 
states appear. We discuss possible implications for cuprate superconducting and related materials. 
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Mesoscopic and inhomogeneous ordering of charges 
in diverse materials such as high-Tc superconductors 
[1], colossal magnetoresistant manganites [2,3], and two- 
dimensional (2D) electron gas systems [4] is a sub- 
ject of intense recent study, as such ordered phases 
may be highly relevant for understanding the properties 
of these materials, including superconductivity mecha- 
nisms, pseudo gaps, transport responses, and magnetic 
responses. Such phases include labyrinths, stripes, and 
clusters, and often consist of intermediately ordered pat- 
tens, lying between completely ordered and completely 
disordered systems [5]. Inhomogeneous charge ordering 
phases can be produced by a competition between repul- 
sive and attractive interactions. In the case of the metal 
oxides, holes with a repulsive Coulomb interaction move 
in an antiferromagnetic background, with the distortion 
of the spins giving rise to a dipolar attraction between 
holes, and allowing the formation of clump, Wigner crys- 
tal, and stripe phases. A variety of other microscopic 
mechanisms that produce a similar coexistence of long- 
range repulsion competing with directional short range 
attraction occur in a wide range of systems, such as de- 
fects or dislocations in elastic media, covalent glasses [6] , 
and systems with a finite density of Jahn- Teller polarons 
[7] . Systems with effective competing repulsive-attractive 
interactions are abundant in soft matter, where bubble 
and stripe phases are observed [8]. A particularly rele- 
vant issue in these inhomogeneous systems is the role of 
temperature. Open questions include how the mesoscop- 
ically ordered patterns fluctuate or melt, and what the 
relevant dynamical modes are. Another important issue 
is what measures best characterize these systems. Tem- 
perature effects are particularly important since dynam- 
ical and morphological changes along filaments couple to 
changes in the material properties of the sample. 

In this work we focus on the thermal properties of a 
quasiclassical model for charge ordering of holes in transi- 
tion metal oxides, in which the particles have a Coulomb 
repulsion and a dipolar attraction. As a function of hole 
density (doping) we observe an extended region that we 
term a "soft phase" comprised of partially ordered fil- 
aments. Ordered clumps form for densities below this 



region, and ordered stripes (Wigner crystal-like phases) 
occur above it. We find that the soft filamentary struc- 
tures persist to high temperatures. Within the soft phase 
region there is a low temperature onset of motion along 
the filaments coupled with fluctuations of the filamen- 
tary structures. The filaments act as a template for cor- 
related percolation of particle motion. When the parti- 
cle positions are averaged over long times, the filaments 
form an ordered checkerboard pattern. We have also 
considered a simpler model of long-range repulsion and 
isotropic short range attraction, and again find evidence 
for a filamentary soft phase, suggesting that the presence 
of soft phases with low temperature correlated percola- 
tion may be a generic feature of systems with competing 
long/short range interactions. However, the filamentary 
soft phase is amplified by the anisotropy (e.g. dipolar in 
this present case), and the dynamics is typically softer 
on filaments than in clusters. 

Our phenomenological model for the holes consists of 
particles with competing long-range and short-range in- 
teractions [10]. Four charge-ordering phases arise de- 
pending on the hole density and the strength of a dipolar 
interaction induced by the holes: a Wigner crystal; a 
diagonal stripe phase; an intermediate geometric phase; 
and at low hole densities and larger dipole interaction 
strengths, a clump phase. The transitions among these 
four phases have been studied in simulations [11,12]. We 
tune the doping level in our model by directly varying 
the hole density in a rectangular computational box of 
size Lx X Ly, with and Ly up to 100 unit cells in 
a Cu02 plane. We initialize the simulation by placing 
the holes at random and assigning each hole a magnetic 
dipole moment of fixed size but random direction. We use 
an efficient Monte Carlo method to find the minimum of 
the total potential in this system [12]. 

We base our model for the interaction between the 
charges on the spin density wave picture of the transi- 
tion metal oxides. Full details of the model can be found 
in [11,12]. The effective interaction between two holes, 
1 and 2, a distance r apart, where r is relaxed to an 
arbitrary (continuous, off-lattice) value, is given by 



1 




600 



400 



200 



(b) 



O ( 



X 200 400 600 800 ^ 



200 



400 600 



400 



300 



200 



100 



(c) 







I o I =1 



o 
o 

o o o o 



o o o 



oo 0° o° o° o" 
? o o o o < 
' o o o o o 

°°i°„°l°o?°o^°^°° , 



400 



300 



200 



100 



(d) 



f^- °o °0 °e 
o o„ o 



X 100 200 300 400 







o o o 

) °n o O 
o o O 



o o " o o o„ o o o 

°„r°„yo °o°°i°°i 



X 100 200 300 400 



FIG. 1. Static positions of holes for different densities, (a) 
Clump phase, n = 0.6, (b) soft phase, n = 1.2, (c) soft phase, 
n = 2.1, (d) anisotropic Wigner crystal phase, n = 2.7. 



V{f) = ^ - Ae-^/" - Bcos(2^ - (/>i - </>2)e-^/^. (1) 



Here, g = 1 is the hole charge, 9 is the angle between 
r and a fixed axis, and (f)i,2 are the angles of the mag- 
netic dipoles relative to the same fixed axis, which we 
assume can take an arbitrary value. A is the strength of 
the short-range anisotropic interaction, and B is that of 
the magnetic dipolar interaction [B w ^/(27r^^)], which 
we assume to be independent variables. Throughout this 
work we take ^ = 0. For the magnetic correlation length, 
we assume the approximate dependence ^ = 3.8/\/n A; 
however, in order to observe a larger number of clusters 
for a fixed number of holes N = 225, we artificially reduce 
the value of ^ to a smaller value with ^' = VO.IS^. 
Test simulations on much larger systems of 2000 particles 
with the full scrc;ening length indicate that the physics is 
qualitatively unchanged by this screening length reduc- 
tion. The distance f is measured in units of the cuprate 
lattice spacing gq = 3.8 A, and the hole density is given 
by n. The sample is periodic in the x — y plane [17]. 

We first consider the static hole configurations, illus- 
trated in Fig. 1, obtained at different hole densities by 
annealing the system from a high temperature. At low 
densities, cross-like clumps form and organize into an or- 
dered lattice structure [Fig. 1(a)]. Above a density of 
n = 0.9 the clumps begin to touch and are replaced by 
partially disordered filamentary patterns of the type illus- 
trated in Fig. 1(b) for n = 1.2 and Fig. 1(c) for n = 2.1. 
These filamentary patterns persist within a glassy win- 
dow that extends up to n = 2.4,when a more ordered 
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FIG. 2. Image of system above the filamentary melting 

transition, with lines indicating the motion of the holes be- 
tween consecutive simulation frames, (a) n — 0.6 in the clump 
phase, where the particles remain predominantly localized, 
(b) n = 1.2 showing the square modulated liquid phase in the 
filamentary regime. 

anisotropic Wigner crystal pattern forms, as shown in 
Fig. 1(d) for n = 2.7. 

The transition from the filamentary pattern phase of 
Fig. l(b,c) to the anisotropic Wigner crystal phase of 
Fig. 1(d) occurs when the interactions between the holes 
become dominated by the Coulomb repulsion. If the sys- 
tem formed a perfect isotropic Wigner crystal, the lattice 
constant of this crystal would be rwc = /VN. We 
can estimate the density at which the transition occurs 
by a balance of forces, q^/f^c « B/S^. This becomes 
Qq/B = rwc, giving n = i/lOB/ao, or n = 2.3 for 
the parameters considered here, in reasonable agreement 
with the transition observed in the simulations. 

We next consider the role of temperature-induced fluc- 
tuations of these patterns. For the clump region seen in 
Fig. 1(a), the superlattice clump structure remains sta- 
ble up to a very high temperature. The superlattice and 
clumps break up simultaneously around T„i — 930 K, 
where the melting temperature is measured by the 
onset of diffusion (see inset to Fig. 4). In general par- 
ticles are confined within the clump up to this tempera- 
ture, although occasionally a particle can escape from the 
edge of clump and move to another clump. Such limited 
motion of particles within the clumps begins at around 
700K. The anisotropic Wigner crystal phases found at 
densities n > 2.4 beyond the filamentary soft phase are 
stable up to about 700K. Above this temperature the 
entire pattern melts rapidly, although we have not deter- 
mined the order of this transition. 

We now consider the fluctuations within the soft phase. 
Well below T^, we find a modulated liquid phase, in 
which the charges are constrained to remain within a 
filamentary square pattern, but are free to move along the 
filaments. This is illustrated in Fig. 2(b) for n = 1.2 in 
the soft phase, showing the holes moving along a square 
pattern. The pattern is oriented with the simulation box 
due to the residual dipoles, which form because we have 
not restricted the spin on each charge to take discrete 
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FIG. 3. Pattern reorientation measure N,n as a function 
of temperature for different hole densities: filled circles, 
n = 0.75, filled squares, n = 1.2, open diamonds, n = 1.65, 
x's, n = 2.25, and open triangles, n = 2.7. 

values but have instead allowed the angles to vary contin- 
uously. The dipolcs align along ±45°, and thus the sys- 
tem is divided into two equal but interconnected basins of 
magnetization. We contrast the filamentary square pat- 
tern observed throughout the soft phase with a similar 
image of the clump phase, shown in Fig. 2(a) for n = 0.6. 
Here the holes remain localized within each clump. 

The modulated liquid phase remains stable up to the 
temperature at which the pattern itself is destroyed. 
We can define a second, lower melting temperature Tg 
at which the onset of the square (checkerboard) fluid 
phase, or modulated liquid, occurs. In this modulated 
liquid phase the particle motion is constrained to follow 
the square template. Although the hole density is not 
high enough to create this square structure in the static 
configuration, when the holes are moving their effective 
density increases, allowing them to form the square state. 
At any moment the square state is not fully formed; in a 
snapshot of the system, the square state contains imper- 
fections and density modulations. On average over time, 
however, the square state is present, and the particles 
move in response to its structure. 

The temperature Tg at which the square modulated 
liquid appears varies with doping. More disordered hole 
configurations (such as n = 1.2) have a lower T, than 
more ordered configurations (such as n = 2.7). To quan- 
tify this, we must detect the transition to the modulated 
liquid. The onset of the modulation occurs when the 
structure of the pattern begins to change globally due 
to local rearrangements of the holes. For example, for 
n = 1.2, one of the loops may open and reconnect into 
a new loop formed out of a previously straight segment. 
We can detect changes in the structure of the pattern by 
monitoring the local orientation (in the range [0, 7r/2)) 
of a Voronoi polygon centered on each particle 
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FIG. 4. Onset temperature Ts of modulated square liquid 
state as a function of hole density n. Inset: Melting temper- 
ature of pattern T„ (open circles), and Tg (black squares). 

at fixed T. Changes in the pattern orientation produce 
increased fluctuations of the local orientation < -\/©J >. 
We find the number of particles A^„, undergoing pattern 
changes by counting the particles that have < -^/Of >> 
0.27r during a fixed number of MC steps. The onset of the 
modulated liquid should occur when enough particles are 
fluctuating on average to allow the excitations to perco- 
late through the system. If the holes are spaced ^' from 
each other along a line, L^/^' = 38 holes fit across the 
sample length, in NS^'/L^ = 5.9 rows of closely packed 
particles. We therefore select = 6 as the filamentary 
transition, since this would allow one fluctuation at each 
crossing point of the square pattern to extend across the 
system. The temperatures Tg obtained with this value 
of Nm agree well with the onset of the modulated liquid 
phase as observed through animations of the simulations. 
We find a clear decrease of Tg in the soft phase when the 
modulated liquid appears, as illustrated in Fig. 3. 

We plot Tg as a function of density in Fig. 4, as de- 
termined from measurements of Nm- We find a dramatic 
decrease in Tg over a soft window from n = 0.9 to n = 2.4. 
This corresponds to the geometrically disordered percola- 
tion phase. Outside of this soft window, particle motion 
does not occur until significantly higher temperatures are 
reached, especially for n < 0.9. Using different values of 
Nm to determine Tg does not alter the range of densities 
at which the soft window appears; it merely shifts the 
curve slightly in temperature. 

In the soft window, we observe correlated percolation. 
The onset of structural percolation coincides with a highly 
directional softness. This is in contrast to the isotropic 
melting that occurs at much higher temperatures. Thus, 
we have a system which contains a rigid template of soft 
filaments. At Tg, a type of dislocation mobility transition 
occurs in which a bundle excitation becomes delocalized. 
In this case, because the effective interaction between 
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particles along the chain is highly nonlinear, propaga- 
tion of a soliton-like excitation can occur, resulting in 
the time-averaged square pattern. 

In the model considered here, the short-range attrac- 
tive interaction between particles has a strong direc- 
tional dependence. We have also studied a system with 
isotropic short range attractive interactions [18], and find 
that in the stripe phase, melting first occurs along the 
length of each stripe. We have measured the melting tem- 
perature T„i of the three phases in the isotropic model 
by computing the particle diffusion at each temperature. 
Particles diffuse freely within the stripes but do not enter 
the regions between the stripes until higher temperatures. 
In this system the stripe structure is not filamentary, so 
only constrained two-dimensional motion of the particles 
within the broad stripe cluster is possible. Tm is lower in 
the stripe phase; however, we do not find as dramatic a 
decrease in due to the quasi-two-dimensional nature 
of the particle motions within the stripe. 

We now discuss some implications of our results for 
experiments on metal-oxide materials. Our results in- 
dicate that the charge-ordered states should persist up 
to very high temperatures. Signatures of disordered fil- 
amentary states occur at much lower temperatures with 
a transition to a square (checkerboard) state at inter- 
mediate temperatures. However, the short /long-range 
interactions only appear upon (polaron-like) localization 
of holes, which onsets at the "pseudogap temperature." 
Above this temperature, a more metallic electronic state 
is expected. The square phase that we observe is a time- 
averaged, not a static, structure, so we predict that it 
would be detected in experiments such as STM which 
resolve time-averaged images, and it may be a good tem- 
plate for the STM images observed in doped cuprates 
[19,20]. Short-time imaging techniques may reveal only 
disordered structures. The low temperature phase at 
T < Ts is not fluctuating strongly, and thus appears 
more disordered than the square phase at T > T, , which 
could lead to an observable disorder-order "transition" 
with increasing T. Noise mcasiirements could also de- 
tect the onset of the square modulated liquid phase: in 
the clump phase there should be little noise, while in 
the modulated liquid phase, large scale fluctuations of 
the patterns should produce increased noise. It is also 
likely that external fields will easily induce currents along 
the filamentary paths in the soft phase. The soft phase 
at Tg < T < Tm also shows similarities with the in- 
homogenous states observed in manganites between the 
true critical temperature Tc and a higher temperature T* 
at which short-range order first appears [3]. 

Some stripe-based theories for superconductivity re- 
quire fluctuating stripes [21]. In our system, we find that 
the soft regions (which appear in an intermediate dop- 
ing region) dynamically fluctuate at low T. Our soft 
phase could then correspond to the fluctuating stripe 
regime. An important feature of the soft phase we ob- 



serve is that the fluctuations are predominantly on perco- 
lating filaments rather than meandering of the filaments 
themselves. The fluctuating checkerboard state may thus 
provide a good theoretical starting point for introducing 
quantum mechanical effects. 

In summary, we have considered a model 2D system 
with competing long-range and directional short-range 
interactions. We flnd a soft fllamentary or flbrillar regime 
that occurs at densities between a low-density clump 
phase and a high-density anisotropic crystalline phase. 
This soft phase consists of partially ordered filamentary 
patterns of charges. At low temperatures, the filamen- 
tary patterns act as a rigid template for a correlated per- 
colation of particles along the pattern: within the soft 
phase, the onset of motion occurs at a far lower tem- 
perature than that for melting the template backbone. 
(Quantum-mechanical mesoscopic tunneling of the tem- 
plate can occur but at much slower rates.) We have also 
found a similar soft phase for other models with compet- 
ing long/short range interactions, indicating that these 
phases may be general features of such systems. 
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